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paper offers five different models in response, followed by development of a useful general expression for 
the optimal control. 


MATHEMATICAL MODELS 


The basic vibration isolation approach that was presented in reference 3 is developed more 
generally below. It assumes a state-space expression for the system equations of motion, having the form 


x = Ax + Bu + f 


(i) 


Single-Mass Models 

Suppose the experiment is modeled as a mass, m, with position, x(t). Let the space station wall 
have position, d(t). If the umbilicals connecting the wall to the experiment can be modeled by a spring 
and a damper of stiffness, k, and damping, c, respectively; the system can be modeled schematically as in 
figure 1. A magnetic actuator applying a control force ai(t) in response to a control current i(t) has been 
included in the figure. 

The system equation of motion is: 


mx + c(x — d) + k(x — d) + ai — 0. 


In state space notation this becomes, 
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where: x^ = x, x 2 = x, and u = i. This is the model which was used in reference 3, i.e., an absolute-state 
model with a position-plus-velocity disturbance . As shown in reference 3, a quadratic cost function of 
j = 1/2 c (x t w i5 + w 3 u a )dt effectively weights low frequency acceleration disturbances by a factor 
proportional to l/u> 4 . 


The system could alternatively be represented by a relative-state model with a velocity disturbance , 
having the form: 
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where Xj and x 2 are now the relative states x-d and x - d, respectively. The same quadratic performance 
index as before could be used. This second model has the advantage of allowing relative states to be 

weighted, but it does not permit direct acceleration weighting, since x is not represented as a state. 


A third model corrects this deficiency by adding u(t) = i(t) as a state x 3 and noting that 
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. States Xj and x 2 are the relative states defined previously. The equation 
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allows experiment acceleration to be weighted directly in the performance index, as shown below. The 
state equations can now be written as 


x = Ax + bu + f 


(6a) 
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If the performance index is 
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where W n is as defined in equation (5), then p j can be used to weight x; W 12 can be used to weight x-d, 
x-d, and i(t); and w 3 can be used to weight the slew rate u = di/dt. 

Equation (7) can be rewritten as 


00 



where 

Wj = PiW n + w 12 
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u 


(8c) 


= as in equation (6c). 
dt 


Once a satisfactory constant feedforward gain has been determined (by the method to be developed 
in this paper) the control could be implemented according to the block diagram of figure (2a) or, 
equivalently, figure (2b). This third model, a relative-state model with direct acceleration weighting and 
an acceleration disturbance , allows for relative-state weighting, direct acceleration weighting, and 
slew-rate weighting. It has the disadvantage, however, that it assumes the system parameters to be 

known well enough for equation (5) to provide a sufficiently accurate estimate of x. Such an assumption, 
of course, may not be the case with umbilicals. 


Two-Mass Models 

Suppose now that the system is modeled using two masses, with and m 2 representing the 
experiment mass and an attached accelerometer mass. Let k x and Cj be the umbilical stiffness and 
damping, and let k 2 and c 2 be the stiffness and damping of the accelerometer. The model is shown 
schematically in figure 3. 

If the system states are defined as follows, 

t 

x 0 = J [ x l( r ) - x 2( r )l dt ( 9a ) 
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x x = Xj — d 


x 2 = x 2 - X 1 


x 3 = x i = *1 - d 
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then the state equations and performance index can be expressed by equations (10a) and (10b), 
respectively, where u(t) = i(t). 
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J = I J (x T WjX + w 3 u 2 )dt 


(10b) 


In this formulation x is approximated closely by the state x 2 , which is available as the accelerometer 
output voltage. This relative-state model with accelerometer feedback and an acceleration disturbance 
adds to the complexity of the problem. However, it allows the experiment acceleration to be represented 
directly as a state. It is necessary only that the accelerometer mass, stiffness, and damping be known 
with reasonable accuracy. 

A system making use of a load cell instead of an accelerometer, as portrayed in figure 4, would lead 
to very similar state equations. Let kj and Cj describe the umbilical characteristics as before. If a load 
cell of stiffness k 2 separates the experiment mass m 2 from a small plate of mass m x attached to the end of 
the magnetic actuator, then the state equations can be expressed by equation (11). 
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Here the states, performance index, and control are defined again by equation (9). This relative-state 
model with load-cell feedback and an acceleration disturbance also allows the experiment acceleration to 
be represented directly as a state (viz., x 2 ). It has the further advantages that k 2 can be known with 
great accuracy and that the load cell output gives an accurate state measurement even with very low 
frequencies. 

All five of the above models have the form 

x — Ax + Bu + f (12a) 

and cost function 



These equations, along with the initial conditions 



and the terminal conditions 
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have already been treated analytically in reference 3. However, an analytical solution is very 
algebraically intensive. For this reason we seek a general matrix solution that is more amenable to 
numerical treatment, to accommodate better the higher order systems. This paper will use a differential 
equations approach to find such a solution. A state transition matrix approach that yields confirmatory 
results will be the subject of a later paper. 


OPTIMAL CONTROL PROBLEM 

The optimal control problem is to find the control u(t) which minimizes the quadratic performance 

index 
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subject to the conditions 
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x(t) — 0 

(14b) 

f(t) - 0 

(14c) 


where Wj and W 3 are the state- and control-weighting matrices, respectively. 


SOLUTION 

The argument of the cost function J from equation (13) is first augmented by the Lagrange 
multiplier A(t) times the state equation of motion equation (12a), to yield the equation 


J = J Hdt 
0 


(15a) 


where the Hamiltonian H is 


H = - (x T W x x + u T W 3 u) + A (x - Ax - Bu - f) 
2 

It is desired to obtain an optimal solution u = u which minimizes T. 


(15b) 
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The first variation of T(x,u,x) is 

00 
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be , , aH , , aH 
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The control u_which satisfies 6T = 0 minimizes T and provides the desired optimal control. Proceeding as 
in reference 3 the problem solution can be summarized as follows: 
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where A_satisfies the system of equations 
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x(0) = x 
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If A can be found in terms of x_and f, equation (16) will provide an expression for the optimal control. 
The solution of the related homogeneous system 
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Assume for simplicity that the eigenvalues of A are all distinct, so that the Jordan canonical form J of A 
is diagonal. Form a matrix X composed of the eigenvectors of A (as its columns) such that 


J = X~*AX = 


A 0 
0 -A 
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has the negative eigenvalues /x* of A as its diagonal elements, arranged such that < H < - < < °- 


Partition X and X as 
homogeneous equation (19) is 
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, respectively. Then the solution to the 
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(23c) 



is to be solved and integrated to find c 1 and c 2 - 
The solution to equation (24b) is 
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Integrating, 


£, = 2, + J •- A ‘X ( ,r 1) M‘ 
= 1 2 + J 


(26a) 

(26b) 


where the 7 ^ and 7 2i (i = l,...,n) are 2 n integration constants and the integrals are indefinite integrals. 
Substituting into equation (24a) yields 
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Applying terminal conditions on x(t) (eq. (18b)) leads to the result that 80 that 

equation (27) simplifies to 
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and 


¥ = e ~ At J e At X^ ^fdt 

If now we solve for A_in terms of x, we obtain the result 

A = (x 2 iX-')x + (x J2 - 


(28c) 


(29a) 


or, equivalently, 


A = (x 21 X n 1 )x + X[ 2 l) V 1 * J e At X^ 1} fdt ( 29b ) 

where X 22 ^ ^ 1 = (X 22 ~ X 21 X n X 12 ) 1 is the inverse of the lower right-hand partition of X 1 , as 
defined before. Applying equation (16) we have the result that the optimal control is 

u* = (w~ 1 B T X 2l X~ 1 jx + (w“ : l B T 34“ 1)_1 e- At J e At X^ 1} f(t)dt) ( 30 ) 

This is a simple form of the result obtained by Salukvadze in 1961 (ref. 1). If the disturbance f is set 
equal to 0, the optimal control reduces to 


H * = Iw.’bV-'J, 


(31) 


By comparison with the solution to the LQR problem (ref. 4), it is apparent that 


-x 21 x-' = p 


(32) 
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where P is the solution to the Algebraic Riccati Equation (A.R.E.) 


PA + A t P - PBW 3 1 B t P + Wj = 0 ( 33 ) 

Making this substitution for — X 21 X u 1 into equation (30), 

u* = -W~ l B T Px + l B T X^ 1)_ V At J e At X^ 1) f(t)dt ( 34 ) 

If the indefinite integral in equation (30) is integrated repeatedly by parts, the control can be expressed 
as 

u* = -W 3 ' 1 B t Px - W3 1 B t X ^“- 1)_1 £ (-A -1 ^ 1 (35) 

r=0 


Assume, conservatively, that no derivatives of f(t) are available. Then equation (35) reduces to 
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where the “(0)” indicates that only the term for r=0 has been retained and the indicates that the 
resultant control only approximates the optimal. It can be shown readily that 
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is the state matrix for the closed-loop system. 
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Substituting into equation (38) 



= - W g 1 B t Px + W s 1 B T A _T Pf 


(41) 


By using equations (17b) and (21) one can show by simple matrix manipulations that 

A = XjjAX ” 1 («) 

Equation (38) can then be expressed alternatively a s 
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or, applying equation (32), as 
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CONTROL EVALUATION 
Physical Realizability of the Control 

The control (35) is physically realizable if the states and sufficient derivatives of f(t) are accessible, 
if the coefficients are real, and if the higher order terms are negligible. 

It can be shown by a state-transition-matrix approach that the coefficients are always real, whether 
or not the eigenvalues are real. Also, for eigenvalues of large enough modulus ( — A 1 ) r+1 rapidly 
converges to the zero matrix with increasing r, so that coefficients of higher order terms rapidly disappear 
as well. Given slowly varying disturbances the derivatives of f(t) will be small, and the convergence of 
the higher order terms will be even more rapid. 


Transfer Function Matrix 


Neglecting higher order terms the transfer function matrix for the closed-loop system is given by 


X(s) = (si - A) -1 (I - W)F(s) 


(45a) 


where 
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Control Stability and Stability Robustness 

Since the control feedback gains are simply LQR feedback gains the closed loop system is stable 
and enjoys the stability robustness characteristics guaranteed by LQR theory, 1 viz., a minimum of 60° 
phase margin, infinite positive gain margin, and 6 dB negative gain margin. 


SAMPLE PROBLEMS 
First Order Problem 

Let a first-order problem be described as follows: 

Find u = u which minimizes J where 


J = 1 J (wjX 2 + w 3 u 2 )dt 
2 0 


for the system 


x = ax+ bu + f 


(47a) 


(47b) 


with initial condition 


x(0) = x 0 


(47c) 


and terminal conditions 


1 Phase margin and gain margin in the multiple-input multiple-output case are measured by diagonal 
perturbation; i.e., the same change in each channel (ref. 5). 
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t— 00 


x(t) = 0 


(47d) 


lim 
t— oo 


f(t) = 0 


(47e) 


The Hamiltonian matrix is 
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Equations (36), (38), and (41) each yield the result 
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( 49 ) 




for which the closed-loop transfer function is given by 
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With LQR F/B alone the control is 
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for which the closed loop transfer function is given by 

X(s) = 1 F(s) (52) 

s — fi 

Note from (48b) that the modulus of fi is always greater than that of a, so that inclusion of the 
feedforward control term always result in better disturbance rejection than that from control with LQR 
F/B alone. 

Comparison of the performance indices and corresponding to and ^ ^u, respectively, is 
very messy algebraically. However it can be shown, for Xq = 0 and arbitrary dwindling function f(t), that 
< J L if a < 0 (i.e., if the system is open-loop stable), and^that y'j < J L if a > 0 if -/i > 4/3 a (i.e., 
if Wj/w 3 > 7a 2 /9b 2 ). For x 0 * 0 and f = 0 the two controls u L and ^°^u are identical, as are their 
respective performance indices. 


Second Order Problem 

Suppose now that the solution method is applied to the space-experiment disturbance-rejection 
problem raised previously in this paper. Consider the first mathematical model of the system, so that the 
problem is as follows: 

Find u = u which minimizes J where 

j = l J (w T WjX + w 3 u 2 )dt ( 53a ) 
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Let the system parameters be m — 100 lbm, k — 0.3 lbf/ft, and c = 0.000622 Ibf-sec/ft, corresponding to 
damping ratio £ = 0.1 percent. 
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Selected results are found to be as previously reported in reference 3, where this particular problem 
was treated in detail. There it was noted that the performance index weights low frequency accelerations 

1/2 

> so that low frequencies are more heavily weighted, as 

u/* a/ 

desired. The result is that for low frequencies the inclusion of feedforward in the control can lead to 
acceleration reductions orders of magnitude below that afforded by simple LQR feedback alone, without 
leading to poor performance at higher frequencies. Different frequency-weighting schemes, of course, 
would lead to different disturbance-rejection characteristics (ref. 6). 


For any given set of weights w 11? w 22 , and w 3 , the closed loop gain for the system with LQR 
feedback alone can be compared with that for the system with LQR feedback plus proportional 
feedforward, at a given frequency. Table I presents the factors by which the DC level (o; = 0) for the 
closed-loop system is reduced by merely adding a proportional feedforward term to standard LQR 
feedback. This represents, in general, the factor by which the closed-loop Bode-o plot for a given set of 
weights is lowered at low frequencies by the addition of proportional feedforward. Note that for all 
combinations of weights considered — each weight being allowed to vary over four orders of 
magnitude — the reduction factor lies between zero and one, varying from 6.00xl0~ 4 to 9.86X10" 1 . The 
optimization procedure accomplishes this task by moving the transfer function zeros to reduce the gains 
in accordance with the frequency-, state-, and control-weightings present in the performance index. The 
poles of the transfer function depend on the feedback gains alone, and do not vary with feedforward 
changes. 
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CONCLUSION 


This paper has developed a general expression for an optimal control in the case of system state 
equations with dwindling forcing (disturbance) terms included, given a quadratic performance index. The 
control, when expanded in series form, has been found to entail constant, real-valued feedback gains 
identical to those determined by the standard LQR approach, along with constant, real-valued 
feedforward gains premultiplying the disturbance terms and their derivatives. 

It has been found that the control offers significant disturbance-rejection improvements over a 
control that uses LQR feedback alone, without sacrificing robustness. In at least the 1st order system, 
sufficient conditions have been presented for the simplest F/B plus F/F control which result in a lower 
performance index than with LQR F/B alone. With large enough closed-loop system eigenvalues and 
slowly varying disturbances, the conclusion was made that only a few feedforward terms are needed to 
approximate closely the actual optimal controls. 

Five mathematical models of a one-dimensional disturbance-rejection problem were suggested, each 
of which is in a form amenable to the optimal control approach presented by this paper. Application was 
made of the optimal control method to one of these models, leading to the same numerical values of gains 
(both F/B and F/F) as found previously using a nonlinear-algebra formulation (ref. 3). 
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TABLE I. - DC GAIN REDUCTION FACTORS WITH THE ADDITION OF PROPORTIONAL 
FEEDFORWARD, FOR A SINGLE-MASS SYSTEM WITH M=1()0 LBM, 

K=0.3 LBF/FT, AND C=().()0()622 LBF-SEC/FT (£=0.1%) 


W^O.01 


W,=0.1 


W 3 =1 


W 3 =10 


W 3 =100 



2 




















Figure 2. — Controller Block diagrams for constant feedforward gain. Both figures (a) and (b) are equivalent 
implementations of the constant feedforward gain controller. 
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Umbilical Experiment 


Accelerometer 



Figure 3. — Physical representations of modeled systems using 
two masses and m 2 , representing the experimental mass 
and an attached accelerometer. 
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